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ABSTRACT 


Using the concepts of fracture mechanics, we develop a theory of the 
earthquake mechanism which includes the phenomenon of subcritical crack growth. 
The theory specifically predicts the following phenomena: slow earthquakes, 

multiple events, delayed multiple events (doublets), postseismic rupture growth and 
afteralip, foreshocks, and aftershocks. The theory also predicts that there must be 
a nucleation stage prior to an eartlKjuake, and suggests a physical mechanism by 
which one earthquake may 'tri^er' another. 

These E»*edictions are obtained by combining two fundamental concepts. The 
first is that 


k = C AttS" 
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and the second, that 





where k is the stress intensity factor, At stress drop, X rupture length, X rupture 

velocity, C a geometrical factor, and V^, and n are material constants. The 

first is a fundamental result of fracture mechanics; the second describes stress 
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c^rosion cracking, a w«U est^ltehcd ph^ical procMS that results in suberitieal 
crack growth. 

We investigate in detail two phenomena of ^>ecial intemt and which are not 
predicted by ordinary fracture mechanics: nucleation and delayed multi^e events, 
ht the first case we find that all eartinjuakes must be prece<ted by quasistatic ^p 
over a portion of their rupture surfaces, but it may be difficult to detect in 
practice. In the second case we studied two pairs of delayed multiple events that 
were separated by the same *barrier* in order to caleidate n. We find tittt the 
stress corr<»ion index, n > 24. 


GLOSSARY OF SYMBOLS 


C ■ geoMtric co&ttaat • 1 for « t«o-diaen«ional crack 

■ 20/7ri/2 £or a circular shaar crack 

k • stress intensity factor 

k|j ■ dynaoie stress intensity factor 

k • static stress intensity factor 

s 

R ■ stress corrosion liait 

o 

■ modulus of cohesion in the presence of corrodent 

R* ■ modulus of cohesion under corrodent-free conditions 
c 

e 

R^ ■ a fixed arbitrary point on the (k-R) curve 

L • barrier width 

” lengths of rupture zones of two events in a pair of delayed 
multiple events 

m ■ frequency of aftershock occurrence 

n ■ stress corrosion index 

t ■ time since occurrence of main shock 

t^ “ delay tim between a pair of multiple events 

t^ “ maximum of rise times of individual events in a sequence of multiple 
events 

t^ * time to failure (instability) or nucleation time of an earthquake 
t' ■ t - At 

u • slip on crack plane 

V. ■ crack growth velocity at the stress corrosion limit 



• t • cMrdituita eorrttpondiag eo 
■ « tutcrial proparty 

% • crack length (for two<-di^n8i(mal cracks) and crack radius (for circular 
shear crack) 

X ■ crack edge velocity 

X - initial X 
o 

X^ ■ X at time t before failure 
8 ■ an elastic wave velocity 

dx m static stress-drop 
At “ short tiaie before instability 




INTRODUCnON 


In recent y^urs the devel<^m«tts of Uneu eUuitks fracture mectanto tetra 
be«i applied to «i io^<»rtant pr^tem in {^ophystest developm«it of a dyttamlc 
moctel of rarthqiMkra [Koetrov, 1966; lUci^ds* 1976; Ai^tremt 1976; Foeaum 
Freund, 1976; Das and Aki, 1977 a, b; Freund, 1979] . The central concept of 
fracture mechanics, which has its roots bi the Griffith «»rgy tMdaime [Lawn aiui 
Wibthaw, 1975] , ^ that for a crMk in an Mastic medium, no ^ro^igation takes 
place until the stress intensity factor at the mmck t4>, k, reaches a raliM a 
property of the medium. K is called the "mo<hilus of ce^si<m" (Kostrov et aL, 

C 

1969] . When k > K^, the Griffith iistat^ty ar^ cmd die crack pro{^;ates 
dynamically with a velocity limited by an elastic wave velocity. 

Most mcides and silicates, howevw, exhibit mwe complicated behavior due to 
environmental effects. For these materials, the crack will {^ocM^te when 
k > K^, wfwre < K^, at a velocity X which is a well <tefii^ functioi of k. 
This propagation is stable aiKi quasistatic smd is referred to as mdxsritical crack 
growth. (We shall use the terms ’Stable", "quasi-static” and "subcritical" to mean 
propagation at velocities much less than the sonic velociti« of tlw medium). Thte 
behavior results from stress induced corrosion at the crack tip, the E^^incipal 
corrodent for the present application being H 2 O. This behavior has been firmly 
estalMished in the laboratory for Mode I (tensile) cracks in a witte variety of 
materials includii^ silicates and silicate glasses [see, e.g., Scholz, 1968a, 1972iq 
Martin, 1972; Wiederhc^ and Bolz, 1970; Atkinsoi, 1979; Lawn and Wilshaw, 1975; 
Knott, 1973] . 

e 

The form of the relationship between k and X does not vary significantly with 
the material; only the parameters in the law vary. As an example we show in 



F^re 1 data on miborltical eraek growth In quartz [aft«* Atkinson, 1979] . The 
empirical relatiowhlp is found to take tte f<^m, 

k ^ ( 1 ) 

or 

X - e*' (2) 

SiiMe n, the stress cwrosion index, is large (12.S in Fig. 1), 1 and 2 are nearly 
indistinguishable. 

SiiKse a stellow focus eartiKjuake is a ^ar crack growing in silicate in an 
aqueous ravironment, eqn. 1 or 2 should be used as a complete de^ription of the 
fracture process. SiiKse an earthquake is a mixed Mode II and m shear crack and 
data are only available for ttra Mode I case, ths involves an assumption: that the 
fwm of the law (txit not necessarily the parameters) does not depend on mode 
(Atkinson's method may actually put tl« crack into mixed Model and m, an 
unst^orted statem^t in Evans [1972] beii^ the only argument to the ccmtrary). 
Later we shall determine n fw an earthquake and show that it is in remark^le 
agreement with Atkins<m's results. Although we cannot prove this assumption, it 
seems entirely reascxiable since when <xie considers tlw physical mechanism of 
stres cwrosion there seems to be iw physicochemical reas<m why the (Mrocess 
slM>uld depend on mode. The Indirect evidence in supp<xrt of this is that rock 
exhibits dilatant creep and static fatigue, both {Mrocesses that result from stress 
cmrrosion, in compre^on and under hi^ confining pressure [Scholz, 1968a; Kranz 
and Scholz, 1977; Kranz, 1980] . We shall also assume that there b a lower limit, 
K , such that when k < K no crack growth occurs. There is only limited data to 



»i^Ft this [Wie<ierhorn and BolZy 1970| Evans, 1972] tnit it well foun<M in 
stress corrosion the<^y. In any ease, our results are not critically dependent upon 
this assumption of a stress corrosion limit, because if it <k)M not e^t, we can 
simply define as a value of k below which tlM eradc velooi^ is vanishii^ly small 
and can be neglected. 

It shmild also be pointed out that in usii^ eqns. <1) or (2) and the fracture 
mechani<» sftproach, we are implicitly assumii^; that the most impc^tant forc^ 
governing the propagation of the rupture are the forces at the crack tip and that 
the friction that acts on the crack behind the crack tip plays no role <m the motion 
of the cradc tip. This is clearly an approximation, but the complete problem 
cannot be handled until the full energy balance, as discussed by Kostrov (1974] 
can be staved. For the pr^ent time, we will have to justify this assumption with 
the success we have, with our present approach, in predicting the observations. 

In this paper, then, we consider an earthquake as a shear crack that 
propagates according to a law given by (1) up to the time when the propagation 
becomes dynamic. This theory predicts a variety of phenomena, all of which have 
been observed for earthquakes. The phenomena which arise quite naturally in the 
theory are slow earthquakes, multiple events, delayed multiple events (dmiblets), 
aftershocks, foreshocks and postseismic rupture extension and afterslip. The model 
also contains specific predictions about the earthquake nucleation process. All of 
these phenomena we now see simply as different facets of the same phenomenon. 

In order to underscore the underlyir^ simplicity of this model, we will first 
qualitatively describe how the various phenomena noted above arise, and then 
follow that with a quantitative model which we will ai^ly to several examples. 


niPLICATIONS OF THE THEORY 

For a crack in an indnito, homogontm^* tUotio me^iun, tlw rtUU^nl^ 
bttween oraek X, and stress Intei^ty faet<«, k, at tl» &Biok tip U 

k « <3) 

regardlen of mode. If k > we immediately see, combinii^ (3) with (1) that the 
crack will accelerate to a catistrc^e (this is a sti'oi^r catastrophe tluui a 
Malthtfflian one since the latter is a simple ejqwnential), the relevant equatioia f<^ 
which are given in a later section. This catastri^he occurs when the mrack 
becomes critical and propi^tes dynamically. We shall show later that this occurs 
at a well defined value of k, which we will call K^. Note that is differ«at from 
K*^, the value of k at which insUd>ility occurs in the absence of stress cwrosion. 
K* is a material property that for these materials can be measured only in a 

corrodent'free envirnunent, e.g., a high vacuum. K , on the other hand, is not an 

« 

independent material prt^rty, it depends on n. We do not assume a value of K 

in our prdJlem, we calculate a value of K^. If K*^, then we need not 

« * 

consider K atalL IfK < K then the instability wUl occur earlier than we 
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calculate, but otherwise our results will be unchained. It is most likely, in fact, 
that K. - K* , since as the <»ack-edge velocity approach^ sonic velocities, the 
crack is pr<^>agating too fast fw the corrmient to diffuse to the cradc-edge and a 
vacuum exists at the tip of the crack [Wiederhorn, 1967]. It may not seem 
obvious at first that this will occur for a shear crack, but for a t<^>ographically 
rm^h surface it is unlikely that shear can take place without some dilation. The 
velocity, at is also important and will be discussed in a later section. 


Hw stTMS inttmlty Uot» in (3) is simply tlw i^Utd stress st infinity (or in 
the ease of « orsde with Meti<m, the strra ^op) multiplied by some geometrlesl 
fsetofs. On the oth^ Imn^ tiM frMture «riterimi (1) <m the mstmrial 

Ipi^rtiM siui n* Sfaiee a fault is not a homoi^ema hotii the i^^^Ued 

stress and hence k, ami the roat^al (wop^tiM 8^ airi n wfll be ftatotions 
Dwdtlcm «i the faidt idane. SiiMse a fault plane oM^stt ot two «vfM« in ooitMt, 
and the t<^<^aphy of wrfaoes is known to be Browidan tteyies ai^ Tt^a^ 197S] 
or fractally Brownian [ Mandelbrot, 1977] , we should expect that the heterogeneity 
of aptdied strnses and material pre^rties Is at least as random as the tc^x^^hy 
and that this heterogeneity exist at all Males. If we com{»re a fault pUuM with a 
mathematically flat ^ane, it is <me of the prop^tin of Brownian surfaces tliat the 
standard deviation from a flat plane increases as x^^^, where x is a scale lei^th of 
the section of the Brownian surface that is samj^ed. We sIxHild therefore emwidM 
K, n, and as random variid>les <m tiie fault plane and the rupture proc^ a 
stochastic growth procen [see, also, the disciwsion in Sciudz, 1968b]. Itwn ttw 
sampling dimension x becomes the earthquake radius <x lei^^ (Bmeiuiitm, and the 
cotwequence of the fault beii^ a Brownian surface is that the U^er the earthquake 
becomes, the larger will be the wavelei^^ and inteiwity of heterogeneities that it 
encounters. The only reason smaller earthquakes a^^ar simpler than larger mies is 
due to the fact that we observe earthquakes with tend-Umited instruments and 
becaiwe the higher frequency waves radated by the smaUw heterogeneities are 
more strmigly attenuated in (wopagating from the smirce to tlw iiwtrument. 

For the purpose of the remainder of this diMussion, however, we will comidw 
only gross heterogeneitiM which we will call terriers, after the usa^ of Das and 
Aki [ 1977b] . These are r^ons <»i the fault plane that are particularly resistant 
to slip either because of low applied stfMS a* exceptionally high strei^th 
properties. The term barrier, then, refera to a hetercyeneity of sufficient size that 
its effects <m the rupture pr<^>agation can be d»erved iiwtrumentally. 
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Wt art now rtaify to ditouaa tlit broa<k»r li^to«tl<mf of tte moO^ fVt arffl 
<ki io by (fiaounim aa Mtnariot tfw various poni^ ^miottima ttwt can r^adt 
to spatial ami ttm^^ v^ntlmis in k, and a mi tht fiudt. Th^ Mmitfloa sm 
U liatratfd in Figurts 2 and 3. 

Hu<dtati<m. Hit loadiqi that ia implted in tlw tUotie r^^aid ttw^ 

19103 ia a ttctonic prootn emaristiiv ci a stMdy inmtaat in tht a^td at 
a vary alow ntt, «ieh stitn bdi^; rtitaatd by tht ttftiiqittkt. Hiis ttetonio 
loading la'ootaa ia tqidimimit to a attady InertaM In k f m a pottntial^r growiagf 
crack. It is laually thm^t that no motion occurs until k » whtn tht 
earthquake initiates. With a fracture procMS such as described by (1), however, 
this is impossible. Instead, propagation of the crack begins when k * K^, and it 
quasistatically accelerates up to sonic velocity at k » (Figure 2). It is thus 
fundamental to this model that an earthquake be preceded by some precursory slip. 
We call this stage the nucleation {^laat. The size of the nucleation region and the 
time scale of the process depend only on n and and their spatial distributions on 
the fault. Because of its importance to earthquake prediction, we will discuss this 
quantitatively later. See Smith and Wyss [ 1968] , Sacks [ 1978] , Sadcs et aL 
[ 1980, 1981] , and Kanamori and Cipar [ 1974] for possible observational examples. 

Slow earthquakes. Although the catastrophe implicit in (3) and (1) is very 
strong, there is a finite probability, because of the heterogeneity of the fault, that 
the rupture will propagate into regions in which k < during the nucleation 
phase and stop. What results is a slow earthquake ( Kanamori, 1972; Kanamori and 
Stewart, 1979; Sacks et aL, 1978] . We expect, of course, that this is an uncommon 
phenomenon but worthy of study because it yields a minimum estimate of the 
moment in the nueleatimi phase. 

Foreshocks. During the nucleation phase, we can also expect it to be likely 
that k > K^ fw small regions of the fault within the nucleation region. These 


rtgions will grow dynamically, only to stop by running into adjoining regions where 
k < Kg. When this happens, a foreshock results. Because of the accelerating 
nature of the nuoleatlon process, the pr<Aability for foreshocks to occur increases 
very rapidly as the time approaches the time of the main shock. Note also that 
foreshocks are not an intrinsic part of the nucleation process so that they are not 
required to fit any regular pattern n« are all earthquakes required to have fore- 
shocks. All three of these predicted properties of foreshocks mentioned above are 
confirmed in the observations. This sumests that the study of foreshocks may 
provide inform ati<m concerning the spatial and temporal development of the 
nucleation phase. 

Stooping . In nucleation, only two possibilities can occur; a slow earthquake 
OT a 'normal' earthquake (Figure 2). In the stopping process, however, mwe 
possibilities can occur (Figure 3). When a rupture is propagating dynamically, the 
stress-interaity factor at its tip is a dynamic one, kjj. (All k’s up to now were 
static stress-intensity factors.) The rupture will stop propagating at a point on its 
perimeter when k^ < K^, but slip will continue within the rupture perimeter as 
the displacement field tends to static equilibrium. After static equilibrium has 
been reached there will be a static stress-intensity factor k^ at the crack tip, 
where k^ > k^^ [see Achenbach, 1973, eqn. 5.7] . It is this extra complication that 
produces the additional phenomena. 

For simplicity we will assume that the upture stops at a barrier where 
kp < Kg (a barrier is indicated by a sawtooth in Figxire 3). What happens next 
depends on the properties of the barrier. If: 

k < K . The trivial case results. The earthquake simply stops and the 

5 w 

barrier becomes the end of the rupture (Figure 3a). 

kj 2 Kc > kjj, the barrier is breached before slip stops within the peri- 
meter and a multiple event occurs, with a delay time t^ < tj^, the rise time 
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(Figurt 3b). By rise time we mean the average tUpping time within the rupture. 
This to t)w type oi muitiptot event ^nt dn«lb^ by ^yu and %taie C IBdS] ami 
moftoUed by Das and Aid [ 1977b) . It to worth remarkh^ that because Iwtero- 
gmelty exists at all MalM, eartlxpiakM are their my nature ii^nitely 
multiple. It is simply that our data allow iw to desoribe only the gross 
heterogeneities. 

< kj < K^. This is an interesting case because, unlike those discussed 
above, it cannot be explained on the basis of ordinary fracture mechanics. In this 
case the rupture subcritically propagates through the barrier, going critical when it 
breaches the barrier after a time delay » tj^. This results tai an earthquake 
occurring just adjacent to a previous earthquake and a short time after it. We call 
this process a delayed multiple event (Figure 3c). That this type of event occurs 
has been very well documented for the Naidcai trough of western Japan by Ando 
[1975] . The most prominent events discussed there are the Ansel I and n (1854) 
events, adjacent earthquakes with a time delay of 32 hours, and the Tonankai 
(1944) and Nankaido (1946) events, separated by two years. We shall model these 
two pairs below, and calculate n and from then. This type of multiple event is 
common in some regions [Sykes, 1971} McCann, 1990} Lay and Kanamori, 1980] . 

There are two differences between these two types of multiple events. The 
first to the delay time, the other to more nibUe. Dem>te by the lei^fth of the 
rupture just as the barrier is encountered, and distance the rupture 

propagates after breaching the barrier. In the ordinary multiple event, the region 
has not come to static equilibrium when the barrier breaks, so that the region in 
continues to slip as the rupture {>**opagate8 to its final dimension «i * ‘r and 
the source parameters are those of a single earthquake of rupture dimension 
t ^ in the case of the delayed multiple event, however, the region 

comes to static equilibrium before the barrier is breached and the static frictional 
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•^sgth it rt-ftttbUMM<^ i neft at ti with tlm* d ttatio e«itMt {Sel^ and 
Eiq^Mtr, 1976| DitUrioh, 1971] . In thii em, whtn tht btiritr It b r ttc h ad , only 
tilt fMw rtgi<» fi4M, and tht rwNiIt It two earth^MkM with toiwM pararattm 
i^r^^tt for n^ttvt dim«^« t|tnd tj, 

It it of ootmt ml ttet tht tavritr bt br tMh t d in tIA nm. If tht 

w«w to, ii^Mte in to^tigth with ditttnet fM^ thin /% tht n^tiwt 
iiM^ grmv into tht Nmiiw for tMit and th«i ttop, whm k < m hi 

F^pvt 3d. If ttUt oeourt wt wiU obt^w Dottttiiniie rupturt growth. A numbw of 
OMt of ttdt hatt btM dooufnmtt<^ tht moat r^^Uimt bttaig tlM 1946 Nankai^ 
tarthquakt t An^ 167S] , in wM^ half Uit r^tiwt arta falltd (piaaittatksaUy. (It 
it intMtttii^ that tht rtgi«i that n^turtd qtu^tati^y alw h^ iw af ttrthocda, 
Oiv mo(M would ^w^it ftwtr afttr^radcs, Iwt not a eomptott of 

aftfftfMdct. Tht ^ta ^louM bt rt-«cainintd to ^toifioally addkta thit qMttion.) 
If tht M^ginal n^twt arM alipt ^ji^tatteally aa tht n^turt groira <^Mai- 
atati^y, wt would obt^vt aft«^to. Uda apptan to bt a faMy e<minon 
^t of ttooMlary iniportanet, ^net it n<wmal2y OMitaita iMt than 5% 
of tht moment of the main shock [Soholz, 1972b] . bidttd, the Ntnkaido eatt it 
proN^y an txetptional ont, bt^tm n^turt Mnta of large tarthqt^kM art 
eoromonly oba^td to Marly id>ut t^ri^ 1971] . Tht« pottMiamio rt^urt growth 
tmially ^ptart not to extti^ tlM n^ttwt more tlm a troall ptretnttgt of ita 
original length, though txetpti«ial omm may oeoiw and tinw bt of intertat to thott 
who ttudy Miamie g^a. 

Aftwrtwdct. ^ of a hettrogmeotw fault <hring an ttfthquakt will Mly mi 
ai^igrt tend to tht atatk aUp ^atrttiutton n^tottd for Um» hMnogtn«o<it eatt. On 
a loMl tealt, tmall patehea may sl4> 1 m (or more) tlmn mnroundUig rcgi<mt and 
tlwa bt dynamteally loaded, rathiff tim imloaded, diwing tlM tarthqMkt. Note 
that the barrier of Dai ami Aid [ 1977b] it an extreme eatt of tlwtt patehM, i.e., 

. 
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QOe that does not slip at aU. Since the etftlKiiiakt oeeinn <!^nainieaUy, this iMdii^ 
very rapid, and the k for th^ reg^oia Udce uiy value k < K^. Thus the 
initiai conditions are set at the time of the main shock. Any patch which 
Kq < k < will grow quasistatioally to failure with a time <teiay that can be 
calculated from (1) (see Empire 3e). Thus aftershocks are j^edicted by the model, 
and shCHild have the followii^ characteristics: 

a. Since large earthquakes can be expected to be heterogene<HJS, the 
occurrence of aftershocks, unlike f<Mreshocks, should be nearly ubiquitous. 

b. Aftershocks sh<xild be distributed all over the plane of rupture, not 
nece^arily unif<x‘mly, and it is likely that a concentration of them will occur near 
the ends of the rupture, where the large scale stress c<mcentrati<m exists. 

c. If k is distributed randomly between and for the population of 
patches, then the frequency of failure of these regions, m, will be a function of 
time, t, after the main shock and will take the form. 


m - I 
m * — 


The derivation of (4) is given in Scholz [ 1968b] and follows from a static fatigue 
law of the type found by Mould and Southwick [ 1959] . This law was shown to be 
derivable from eqns. 1 or 2, by Wiederhorn and Bolz [ 1970] . 

d. Aftershocks are a second order effect relative to the main shock since 
they result only from deviations from the mean, so the sum of their moments 
should be only a small fraction of the moment of the main shock. This is one way 
to distinguish them from delayed multiple events. 

One other property of aftershocks which is not directly implied by the model 
but seems likely is that aftershocks within the perimeter of the rupture zone will 
statistically tend to occur in isolated patches and hence these aftershocks will not 
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ttnd to havt aftershock sequences of their own. Aftershocks that occur or tlw 
perimeter, and which theref<»« extend the perimeter, may, the other hand, 
produce after^ociu. This effect wiu ^>served by Pa^ 1 1968] for the 1964 A^dca 
earthquake. 

Hie propertiM of aftmhocks are sufficiently well known that it is not 
necessary to cite particular examples to state that the observations bear out the 
above ii^edictioia of the model. The prediction of tl» Omwi Law (eqn. 4) is 
particularly important. In order to see how coiuistently aftershock sequences obey 
this law and the other sd)ove predictioiu, one should cduult tte compendia of Utsu 
11989, 1970, 1971, 1972]. 

Deep earthquakes. In the above discussioia we have tacitly been ccmcerned 
only with tectonic earthquakes of shallow focus, as defined in the usual way. Deep 
focus earthquakes, on the other hand, do not conform to some of the predictitxis of 
the modeL Most prominently, they do not have aftershock sequences. This may 
result either because the corrodent responsible for this behavior, free H^O, is not 
present at those depths or simply because the mechanism of deep earthquakes is 
not the type of rupture process that we are discussing. 

In the above we have taken a time dependent fracture criterion, which is well 
established experimentally, and applied it to tne earthquake process. Since we 
know that always, kg > kp always and kjj < when the crack stops 

pr<^agating dynamically, every possible relative condition between the k's and the 
K's has been considered in the ^ve discussion. We f<Hmd that in so doing, the 
theory predicts all of the many facets of ruptwe in the earth, many of which had 
no pricHT explanation and further, that no phenomena have been observed that are 
not predicted. We thus now have a physical basis for understanding these 
phenomena. We can see why foreshock sequences may be variable, whereas 
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aftershock sequtncM are very rq^ular. We have a ph^^ieai expAutttimi t» the 
Omori Law, and we can see that since delayed multiple ewnts are a fecial ctM 
and slow ear^uakes, an even mwe q>ecial case, that these crtwrom^tt should be 
uncommon. The success of ttw theory in predictii^ the otaervationa demonstrates 
that rupture in the earth a la# similar to that observed in the labOTatoi^ to 
result from str^ corrosion. It not, of courae, prove that the causative 
mechanism is nec^svily the same. What must be emphasized is that suberitical 
crack growth, governed by a rate equation of a type similar to eqn. (1) or (2), must 
play an important role in the earthqui^e mechanism, altlH)ugh the underlyii^ 
mechanismis) cannot with any surety be i<tentified, and in principle may be 
unidentifiable. 


pppi 


t 
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»LIIlMlkTICAL 

We dei^<^ here a my simple the<mti^ ap^Mch to <»leulatt pr^^lM 
strea- 0 (»rodon ctr«»I^S. Tbe stress-4ntei^ty factor k o ■ W 

where o > stress at the ^ok t^ somS r » <fistoiee from mratok t^> for a two> 
dm«i^<mal plane endc (<rf any mode) aikl for a drmilar plam shm eiadc, ki mi 
infinite homogeneous medium which is linearly elastic everywhere off the <^ok 
plaiw, is given by 

k = C At /X" (5) 

where C 's & geometric ctmstant. At is the static stress drop and X is the crack 
l«^;th for a two-<fimensional crack aiul ^ the radius for a circular crack. C is 
equal to 1 for a two^mensional shear <x t«isile crack [Lawn mid Wilstaw, 1975] . 
For a circular stwar crack, there are two mo(tes ( viz. Modre n and m) of 
prope^ation at a point along the crack e<%e, each with its own stress-intmisity 
factw. Takii^ k to be the square-root of the sum of the squares of these two 
stress-intensity factors, C is giv«i by 20/7 ti/F" [SU i, 1973] . At refers to a tmisile 
or a shear stress component (pending on the mo<te of crack pr<^>e^ati(»i. 

From experimental results on stre»-c(XT(»i<xi crackii^ [Atkins<m, 1979] we 
find 


or 


0 


(6a) 
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k«W^(X)^^” C6b) 

where and are the valuM of k and X at tto initiati(»i of quaaistatie 
growth at time t * 0, and n is a dimensionl^ quantity called the ^^ear‘C<^rosion' 
index. was diseiamd in detail in tiw introdueti<m. The time t a 0 may aetudly 
refer to the time at which the orack velocity is no l<x^^ tnmishii^ small. 
and n are the material properti^ of the medium. Combinii^ (S) and (6a), 


X * V. 


C 




If 


is ind^ndent of time, then X 



(7) 


Integrating, 


X = 



n-2 


V.t 


X 


n/2 


2 

2-n 


( 8 ) 


The free parameters of this equation are and n. The value of n (n is always 

> 2) at this point is simply taken from Atkinson [1979] but later we shall 
determine it. Note that Figure 1 shows that n is independent of temperature and 
humidity but (K^, are not, so that in the earth we expect to find n to be close 
to Atkinson's result but (K^ to be very different, 

Usii^ (8), we can find the time to failure (defined as the time from which the 
crack can be detected to start growing subcritically to when it reaches instability) 
simply by settii^ the quantity in square brackets equal to zero, since this is where 
k goes to infinity, since n >2. This is the 'nucleation time' of an earthquake. Time 
to failure is then given by 


and th# orad( size at time 4t bef<^ failure oan be determined by stdMtitutii% 
t* a tf - At into (8) The vetocity at time At bef^ faUiffe can then be foimd 
A*om (7) by substitutii^ X^^ f«* X. Note tiiat as loi^ as At is in<^>endent of time» 
it (toes not enter into equatioiui (7), (8), and (9). 

At this point let va point out an a<hmntage of toing (9) to determine time to 
failure t^ over the u»ial way of findii^ t^ [Evans, 1972] . t| here depends only on 
the initial conditions and n and on the filial ctmditions. If we now s^ume a 
value for At , we can also calculate K just before failtff e. By fidlowii^ tto above 

C 

metho<j, we have forced the cradc to grow in accordance with a given k > ^ 
relationship (namely, equati<m (6)) until it reaches instability. The e^ve method 
combines the two re-inforcii^ effects of the str^s-corrosion imtability with the 
geometric instability. Note that if the stress-corr(»i<m index n and the point (K^, 
V^) of the (k-X) curve at which the instability occurs are known, the subcritical 
rupture process cannot be determined, but if n and (K^ V^) are known, the total 
rupture process to instability can be completely determined. In this sense, is a 
more fundamental property of the material than K for cases when stress-corrosion 
cracking occurs. 

The formulation developed above is applicable to the case when the crack 
does reach instability. If the crack does not reach instability but propagates 
subcritically through a region of length L, say, then we need a minor modification 
of the method described above. As an example, let us consider the case of a pair 
of delayed multiple events, with the rupture length of the first rupture being X^ 
and with the barrier length between the two rupture zones being the length L. 
Then, combining equations (5) and (6b), we get the delay time between the multiple 


events to be 





W 


*1*1^ <el7> 


a-n 


- (X, + u 




to 

ao) 


The unknown parameters in this equation u>e n, &t , and L. If we have at least 
two se^ of delayed multiple events aeron the same hairier (so that Ly C, &t • 
n are the same) hut the ruptwe lei^in of the first event of each set are different 
and the delay times fw each pair are different, then we shall ffet a seetmd equation 
Ihce (10) with a different values of t^^, say t^ and a different say Dividii^ 
the two equati<ms, we get 


2-n 2-n 

h , X, 2 - (X^ * m 

t- ” 2-n 2-n 

^2 2 
O O 


(U) 


frmn which n can he (Stained provided we assume a value tot L. Once n is ftnind, 
if we assume a At and take to be the strras-intensity factor due to a rupture of 
lei^h X^ we can also determine from (10). In the above we used equation (1) 
together with (5) to derive the mathematical formulation of the problem. We also 
derived similar relations using (2) and (5) and our r^ults of the later sections were 
found to be virtually the same, and so this case is not separately discussed. 

We point out here that this method is valid only up to the point when the 
crack-e<%e velocity approaches the sonic wavespeeds of the medium. Once the 
crack propagates with velocities comparable to sonic, the problem becomes a 
dynamic problem and has to be treated as such. 

In the next section, we shall use the method developed here to model in detail 
some of the scenari<» described earlier. 




1-^ ' - - - .^ . . -^-^.„i«-- ' TniiimWiMnlt-';-*------ 
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THEORETICAL MODELLING OF NUCLEATION PHASE 


In <Mir ealeuUtioni \mixt$ the mettu^ of the prevloua MOtl«it we shall asMme 
that the strength of the material aloi^ the zmie throi^h which the rupture aerially 
^oeeeds to be omstant, Le., and n are cmwtants. Clearly, this will not be true 
in reality but since we do not know the details of this variation of strength, we do 
not think it meaningful at this time to compUcate the model by introdiMsii^ thMe 
variations. Our method, howewr, can be sailed with minor modificatioro to the 
case of variable and n. In that case the fault pUme can be divitted into 
segments of constant and n (concentric on« t<x circular cracks) and our 
method can be applied individually to each section until ii^tability is reached. It 
will generally be intuitively very clear what wcHild happen if were lar^r than 
the one chosen for the calculations. 

We model the nucleation phase of the earthquake as a subcritical exteraitm of 
a circular plane shear crack. The eases we consider are shown in Table 1. X^, 
and n are the input parameters. The values of cik»en were made consistent 
with the given n's, a strws <fr<^ At of 100 bars, and a point (Kj^, Vj^) wi the (k-X) 
curve given by * 10^*^® bar /cm; Vj^ = 0.1 x 10 ® cm/sec for cases (i) and 
(ii), by = 10^*®® bar /cm; Vj^ - 0.1 x 10“® cm/sec for cases (iii) and (iv); 
and by = lo®’®®® bar /cnT, = o.l x 10“® cm/sec for case (v). (How- 
ever, we could also have chosen arbitrarily.) 

The values of U MNm*®^^ = 10^ bar /cm] for cases (i) - (iv) are 
within the rai^e shown in our Figure 1. The different eases represent different 
materials and the time to failure will tell us how stroi^ they are. Using 
equations (8) and (9), the time to failure and the crack-radius X^ just before failure 
can be (Stained. Let us take X^ to be the radius 1 second before failure. The 
results are shown in the last two columns of Table 1. 
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To stu<ty how the eradc ap^oachM instability in e«!h ease, we (dot X n, 
time from equation (8) tw fl^pro^mately the last 100 seomds befwe failure. 1%# 
r^ults are shown in Figur<H 4 and 6. The sh^ of the curves ^(ree with tiMae 
obtained in the laboratory by Wiederiwn { 1967} F^^ve 2] , and, more interMt- 
ii^ly, with the dilatometer records of Sa<dcs et al. [1981] r^ro<teeed in Figure 7. 
The curve of strain as a funoti<m of time jittt prior to an earthqiuke was fmmd by 
these authtm to be exactly of the fm^m shown in Figum 4 and 6. Fmr ease <v), the 
velocity at the last few time steps is clearly comparable to s(mie velocitin of the 
medium so the values of t^ and in this case are not exact. The crack radius 
increases very stowly until the last few time stec» when it increa^ very re^idly. 
It is only during this last phaoe that any preeurs<X 7 strain chaise may be lai^e 
enoi^h to be detectable. However, other precursory phenomena may result 
indirectly from the quasistatic rupture growth diffii^; the nudeation phase. 

In Figures 5 and 8 we plot the k*^ curves alof^ which the crack extends. We 
plot the point and the last three points at the last three seconds 

prior to failure on each curve. If is the stress-intensity fact<» one second prior 

to failure, is given by the topmost point on each (k-X) curve. We have assumed 

* 

K- < K_ in our calculations here. If K. > K„*, the instability will occur earlier 
than we calculated, and X^ will be even smaller and the strain-chaise less detect- 
able. 

Thus our model of the nudeation phase of an earthquake implies that there 
must always be precursory slip and explains why it is seldom detected. We have 
assumed that only one crack is involved in the nudeation process. In reality there 
may be many small «?racks growing subcritically and if enough cracks are involved 
in this process, the resulting strain chaisa may be large enough to be measurable. 

It is dear from Table I and Figures 4 and 6 that the size of the region over 
which the nudeation process occurs prior to the instability as well as the duration 
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of tht nueleation prootss art highly dtptndant on tht valuti of tht variaUta 
anumad f«r tht calculation. Therefore we have no i<tea at ^Ment as to either the 
spatial or temporal scale of the ph«iomen<mf but can oidy make predetitm w to 
its form. The scale can only be determined from observations such as those shown 
in Figure 7. 


THSOR^CAL MODSLUHQ OF DELAF^ MULTIPLB EVENTS 

To s«»oti for this Qrpt ph tnom tnon wt mtot look for two iHtt 

iicwfy i^kstnt wesm thot oeemtd with a tirat <Sefa^ rottoh U»9» than 

the riM Uint of an li»SviAial o^mit but lauoh slwr^ ttan tht ramrrwiot Uaa, 
i. 0 ., 10^ SCO « fj) « Fwfthtrmort, ttw steomi tftnt mu^ bt 

Initiattd in tiw r«^(»i a<|f«^t to tlw n^tin^ of thr fint o^nt, 
pr^agato away hrom it. Hw see«^ en^t of tlw pair is 'tr^^^dF ^ Om flrtt and 
om tlw^ a ^^eal mt<^uuirisin by wM<^ mm aartiapHri^ say * » 4E i w* 

aiK>tisr. 

The Kli PeniiBiula Iwriw. Using an historical record that datM from 
684 A.D., Ando 1 1975] has shown that this phenomena) has repeatedly occwred in 
large earthquakes along the Nankai trough in southwest Japan. He found that the 
plate boundary could be divided into four segments, A,B,C,D, that either rupture 
singly, in adjacent pairs, or all together in a sii^^le great earthquake (Figwe 9). We 
show these regions in Figtire 10 and note that Ando's observaticm indicates that the 
barriers between these regions are persistent, identifiable features. This is 
particularly true of the barrier of the tip of the Kii Peninsula, between regions C 
and B. Since 684 A.D., this bnrrier has been one end of the rupture zone of nine 
l^at earthquake. 

Ando also showed that the Kii Penirwula barrier he ttq>tured with a time 
delay tom time during the historic record, pro<hicing four pairs of delayed 
multiple events. The most reliable record date from 1707, when all four regioe 
ruptured in a single event (Figure 9), thus resetting the initial conditions in all four 
regions at the erne time. The next event we in 18S4, when the Ansei 1 event 
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nipturtd DC, fdtowtd 32 houn Ittw by tht n^turt of BA In tht AiHd Q ovtnt. 
Th# ntxt fvont wu tht Tonankti tvtnt of 1944, whieh rupturtd C, fdlowtd two 
ytm tht Nai9aii<h» irt^ nv«ia>t41 B c^mhwfty and A 

quaiiftatioally. Ando waa abto to show oltirly that tht Naitt^to Mnthi^pttn 
initiattd iHMur tht BC btHoiduy a^ tht wttt, ai prtdiottd. 

Thttt two Mts of dtlaytd multipto ti^nts (1854 and 1944-44) rupttvtd tht 
tamt barritf quaiiatatieally and this providtf us with a uniqut opportunity tor 
ea l cu l athy n for this barritr. Tht rtquirtd modtl paramtttrs for quantitative 
modtlling of thtst tarthquakt pairs art taken from Ando, and stown btlow usij^ 
our ^vious notation. 

Ansti I-n pair Tonankai - Nankal(k> pair 

t^ ■ 32 nours » ?50 da^ 

X * 210 km X*. * 133 km 

o o 

Let us 9dv9 the problem using the various valuM of L shown in Tablt 2. For each 
case tht calculated n is shown in tht stcond column. 

Wt set that the value of n is virtually independent of L. More importantly, n 
differs from the value found by Atkinson only by a factew of 1 or 2. This is a 
remarkable agreement considering the simplicity of our modtl and tht fact that wt 
model a Mode II crack while Atkins<m's crack prop^ated mainly in Mode I. If wt 
•isama At * 100 bars, and take to be the stress-lntenrity factor due to a 
rupt’flre of length .^30 km, wt can determine from (10). 

For the multiple event to occur, the barritr must have a strength such that 

^o-kg* For the eases studied here, X^: lOO km and tt > 100 bars, so that 

3 1/2 2 

10 bar km . The M>«oific fracture energy is given by G * K /4u [Lawn and 

Wflshaw, 1975] , so that G ’ 10^^ erg/cra^. This value, similar to that obtained by 
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Aid [1979] ttit Fwt to ftrf oi 

w&nt tai U)t vidid^ 9t tht »wak waettm wamgf miMh t^wr 

ttwi thftt of ainr ioiowii mturiidt ttwrof ort oart Mould bt takto In ttMrprt^ Iti 
mtanii^. It itoiiit rMional^ t/t> n|^t ttet tl^ to ^l^ttd in tetlatUo 
dtfonnntlMi within a«nt v<dinnt ROTous^tag tht mmek t^ tlte to a Mtiwai 
txtanfion ^ tht 'and aont* eonetpt of Bartnblatt [1999] dtoeuntd ^ Aid [1979] 
in Ml trtatmtnt of tlw tub|tot. A prowiniotd o1m«|^ in tte dip of tM ^d^f 
Mnt uMoiattd with tlw Haidcai trou^ oeotav b«wati> tiw EU Ptnhwda, Imiet It 
to Uktly that tMi tarriw to oauMd by a jog In tht n^tun ^aiw (K. Mogi, ptnoml 
Mmmimieation, 1980). Wt RivMni ttat tht «mgy rtquinKl to atop tht n^ttwt to 
<iiripattd in t)M tMunt bthlml tM jog. 

Ow sitimatt ^ for tht EU P«dmda Mnriw to M^al orctora oi 
mafpUbidt higiw ttwn that muntd in oiw dtoc u Mton oi tht lu^atMi prot^m. 
Wt would U^tfwt not tjq^tot m^tati«i to o<»tff witMn R»h an anmna^ia 
rtgi<nt, at txetnii^ dmt (or strm) wouk) bt rtqtdrtd Variation in E^ or E^ by 
ordtn of magnitudt was also found by Aid [1979] . This suggtsts to us that thtst 
miatimw art <hit to gtomttrieal R>m{dtzitlM of tht fault »int, rathtr than 
changes in tht properties of the fault zone materials which are unlikely to vary so 
widely. It to interesting that, as <Useusstd earUtr with respect to heterogeneity, 
RKdi gemnttrie irrtgulariy to expected to bio’eMt as wMch to exMtly the 
manner In which k scales (at constant This suggesto that the probability for an 
earthquake to stop will be set le independent. 


An Scamole and Coimt^example from the Aleuttaia 
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of I^t ttrtiiqutfcM in tht Atertte k9 

Mt QRf tie« ^ t ttit 1157 Alratkua 

MVtiiqui^ mm fdlowtd 7 ifoift Ut» by tht Rat Mind wl^ initiittd 

It Md c€ tilt IIS7 n^tart »m and pr^i^ttd to tb» mmt 1971] . 

tin tMriff in tfdt oMft in tht mm ot raoit in tm iUrati«t wh«« tht 
n^vt Mhtt lifft tartiiipMc«i ^mteatt^ twi itraif m^tm^on in tht 
bathymt^ [l^«iet, 1977] . 

It <N> flnt b^Mtion, alM i^ir trmo F^urtll ttat tht Not- 

tmbw lOt 1938 tirthqiMlct and tht 1964 AlaMcan t^h^iakt art aimthw txatij^t 
of a Ma^ multi^t tt«it ^dr. Hds to ttw mm^ towtm, bt^at tht 1964 
tv«it initiated in tht northttft part bl ite n^ttat som and propagated to tht 
Muthwtft, Lt.y tommrn tht n^tiat vmo of tiw mtU^ tvtnt. Tiwrtfort tht oloat 
oMioT^iet tai timt oi tiMtt two e«tt%iioia tvmtt to a e<^iei<tonet. 




Wt note mtt additional interMtii^ txam^ i^oni tht Alcutiana. Hm Rat 
Island tarthquakt of 1965 waa follow td a ttqutnet of laigt normal (aultiiif 

I 

tarthquakM in the adfaomt outer wall of tht trtiwh [Staiatort 1968] . This tfpt ot 

ttfthquakt rmdts toom ftexuw ci tht outer wall and rtot rtmilthy f^m 

tubduction [Sykes, 1971]. Tht sudden subduetimi pro^ietd by tht Rat Island 

earthquake would be expected to tn^str* tarthquakta in tht out«r wall bt^tot it | 

would result in a mtddtn instate in the fMcurt and it would also ramilt to a sitodm 

reduetton to tht horizontal oom^taai^ strata in tht outer rte. Thus k would be 

iimttttd f<w normal faults to the outer rtot as a rtmilt of the Rat Island 

tarthquakt, and for any rtgtei to which < k < teiliwt will o^ur after a 

timt delay. Hwst events art U^tfort afttrsiradm of tht Rat Mand tarthquakt, 



28 


but of a differtnt type than disclosed earlier. The Itfgest event in this sequimoe 
was an event of M. » 7.5 which occurred 57 days after the Rat bland eartt^i^ke. 
This is thus a minimum estimate of the nudeati<Mi time t<x tlds earthquake. (It b a 
minimum estimate because our earlier definition of nucleation time wso that it is 
the time to the instability from the initial conifition k ~ K^.) 

Complex Events 

In the above example, the difference between aftershocks and delayed 
multiple events is less distinct than in the usual case. The delayed multiple events 
discussed above are cases in which the two events occurred on the same fault. 
Cfcses of multiple events in which the two events occurred on different faults have 
also been (koerved. Delayed multiple events occur commonly in the Solomim 
Islands, and some of these occur near the sharp comer of the trench near New 
Britain and go around the corner [ Lay and Kanamwi, 1980; McCann, 1980] . This 
the first event has a N-S strike, the second E-W, or vice versa. Lay and Kanamori 
[ 1980] argue that this indicates that the slab is continuous armind the corner, but 
this need not be so. The two Gazli earthquakes of 1976 occurred on two conjugate 
thrust faults with a time delay of 39 days, and certainly in this case an argument 
that they occurred on the same fault would be a forced one [Kristy et aL, 1980] . 
There does not need to be a delay in such complex events, however. The 1927 
Tango earthquake ruptured two orthogonal, cimjugate, strike-^p faults in SW 
Japan. It is not known, however, whether ihis was a ni^mal multiple event in 
which both faults ruptured dynamically or if one fault ruptured dynamically, the 
other quasistatically. We certainly know that it is possible for the first event to 
rupture dynamically and the second quasistatically on different faults, since this 
phenomenon was observed to follow the 1968 Borrego Mountain earthquake in 
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California [Allen et al., 1972] . The earthquake/ which was <hi the San Jacinto 
fault, triggered quasistatic slip <m the neart>y Superstition Hills, Imperial, and San 
Andreas faults. Although we are not presently (urepared to mo<tel thea« complex 
eases quantitatively, we believe they mult from the same phj^ieal mechanUim as 
discussed above. 

The Effect of the Phase of the Loadiiut Cycle 

The loading cycle (<x a given region the time required (<a the tectonic 
process to replenish the stress <fr<^ped in the previoi» earthquake on the same 
secti<m of fault pkne. The period of this cycle is the recurrence time for 
earthquakes in that place. It is clear from the above timt tlw condition fw the 
occurrence of multiple events and delayed multiple events is that the two regitms 
on either side of the barrier must be nearly in phase in their loadii^ cycles. We 
m^iht expect, then, to see adjoining regions going in and out of phase in a cycle 
much longer than the loading cycle. Indeed, two such phase cycles appear in Ando's 
history (Figure 10). We begin in 887 A.D., when an ABCD event occurred, settii^ 
all regions into phase. This was followed by two delayed multiple events: the 1096 
A.D. - 1099 A.D. pair and the 1360 A.D. *1361 A.D. pair. Then the 1498 A.D. CD 
event occurred, with no cwr^ponding BA event, suggesting that BA was too far 
out of phase to be ruptured as a delayed event. The BA region ruptured 107 years 
later in 1605 A.D. The ABCD event of 1707 again put the regions in phase. Next 
follows the Ansei I and n pair, and the Tonankai-Nankai(to pair. The regions seem 
to be getting out of phase again because D did not rupture in the most recent 
sequence. This region (Sun^a Bay and the Tokai district) will either rupture on its 
own, which has not been oteerved to happen befwe, or skip this cycle and rupture 
in the next ABCD or CD event. 




The reason we e:^ct such phase cycles to occur is that two competing 
effects are at work. The first is the rupture proce» itself, which simply by it 
producing rupture on all possible regions in a given event, has a smoothing effect 
on heterogeneity which mults in high phase correlation of at^oining regions. On 
the other hand, spatial variation of the tectonic c»‘oce^ such as variati(m in the 
slip vect<»^ magnitude and direction alot% a plate b<Mn<tary, and heterogeneity of 
the geometry and material properties of the fault ztme have a rmighening effect. 
The first effect produces str<M[ig cliaterii^ in a space**time sense, the sec<md 
randomn^s. What results, and is observed, is weak clustering. This is why space- 
time diagrams, such as those of Kelleher et aL [ 1970] seldom show significant or 
clearly obvious trends. 


31 


CONCLUSIONS 


We have presented a detailed disciosion of subcritieal crack growth and its 
appUcati<m to Utne-depen<tent rupture in the earth. Usu^ very basic ctmcepts of 
fracture mechanics and results (k>tained from lab<^atory ei^riments, we f<Mind 
that our theory predicts various scenarios for time-dependent rupture in the earth. 
We have pr^ented several example to show that all of the (^edicted p.ienomena 
are actually observed in the earth. The theory does not contain any E^edictions 
other than those observed to occur, n<»r do the observations indicate that any 
phenomena occur other than those predicted. Our theory explaii^ why some of 
th^e phenomena are more common than others. Using a simple theoretical 
development, we model two cases in detail, the nucleation stf^e before an 
earthquake and delayed multiple events. For the nucleation problem we show that 
all earthquakes must have precursory slip but the r^ulting strain changes may be 
much too small to be detectable. For delayed multiple events occurring aloi^ the 
Nankai troi^h, we f<xmd two pairs of delayed mutUple events that were separated 
by the same barrier. We uniquely determined an estimate of the str»s corrosion 
index for the barrier between these events. We thus obtained estimates of the 
material prc^rties of barriers in the earth. Our theory also suggests a physical 
mechanism by which an earthquake can 'trigger' another earthquake on an adajcent 
zone of the same fault or on a different fault in its vicinity. 
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1 Sacot^ Bafora 

lndox» 

Radius, 

Valocltjr, 

Tlas to Failure, 
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12.5 

0.1 ea 

0.169 X 10-7 

1125906 sac «. 13.0 days 

1.76 ca 

20.0 

0.1 cm 

0.146 X 10-® 

7593035 sac « 87.9 days 

0.63 ca 

12.5 

0.5 ca 

0.395 X 10‘® 

24094025 sac « 278.9 days 

14.18 ca 

20.0 

0.5 ca 

0.143 X 10"® 

388763365 sac • 12.3 yr 

5.69 ca 

12.5 

1.0 ka 

0.1099 

173382 sac » 2.0 days 

1.4 ka 
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we have eases A, 8, or C. If the barrier increases bi strength with 
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a rcMlt of non-uniform dip on the faidt. 
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shown. The arrows Indicate that the second nature initiated from 
the end where the first rupture tmd been arr«ted by the barrier. 

Fi'^jre 10. Bathymetric map off the smitheast coast of Japan showing the 
segments A,B,C,D aloi% the Nankai trm^h. The terriers between the 
segments may have some weak correlation with the bathymetry. 

Figiffe 11. 'lap of rupture zon« and tlw space time sequence of lar^ earth- 
quake aloi^ the Aleutian are [after Davie et aL, 1980 and Syke et 
aL, 1980]. We make specific reference to the adjacent 1957 and 
1965 events, which are delayed multiple events, and the 1938 and 
1964 pair, which are not. 
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